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I We calculate the ultrasonic attenuation in magnetic fields for superconducting states with line 

■ nodes vertical or horizontal relative to the Ru02 planes. This theory, which is valid for fields near 
^ ' Hc2 and not too low temperatures, takes into account the effects of supercurrent flow and Andreev 
^ , scattering by the Abrikosov vortex lattice. For rotating in-plane field H{Q) the attenuation a(0) 

exhibits variations of fourfold symmetry in the rotation angle Q. In the case of vertical nodes, 
the transverse TlOO sound mode yields the weakest (linear) H and T dependence of a , while the 
longitudinal LlOO mode yields a stronger (quadratic) H and T dependence. This is in strong contrast 
to the case of horizontal line nodes where a is the same for the TlOO and LlOO modes (apart from 
a shift of 7r/4 in field direction) and is roughly a quadratic function of H and T. Thus we conclude 
P3 . that measurements of a in in-plane magnetic fields for different in-plane sound modes may be an 

O ' important tool for probing the nodal structure of the gap in Sr2Ru04. 
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^ ' I. INTRODUCTION 

Q . A number of experiments give evidence that the superconducting state in layered Sr2Ru04 u consists of spin-triplet 
Cooper pairs with broken time-reversal symmetry. rMore recent experiments at low temperatures have established 
power law Tj-dependence of the specific heat (T^), l3 the spin-lattice relaxation rate (T^), S the electronic thermal 

■ conductivity,aiJ the penetration depth,! and the ultrasonic attenuation.LrH These properties are most naturally ex- 
plained in terms of a spin-triplet order parameter d(p) — Az{px + ipy)g{p) where the even-parity function g{p) 
has vertical line nodes (e.g., p^ ~ Py or PxPy ) or horizontal line nodes [e.g., cos{cpz + clq)] on the cylindrical Fermi 

, , ^ surface.Lrilj Since the measured anisotropy of the in-plane thermaLconductivity k[0] for rotating in-plane magnetic 

00 field □ is smaller than the calculated anisotropy for vertical nodes JiJ and since the anisotropy of the inter-plane k[8] 
is insignificant Jj these authors have discarded vertical nodes and suggested instead horizontal nodes in the supercon- 
^""^ . ducting gap. Receittiy we have shown, however, that the amplitudes of k[8] for vertical and horizontal line nodes 
' are about the same Jl3 The small size of the observed amplitude of k[8] is due to the result that the amplitude of the 
variation of k[0] decreases with increasing impurity scattering and temperature. The data for ultrasonic attenuation 
a in Sr2Ru04 measured for different sound directions q and polarizations e in the ab plane are found to be consisteri 
with a vertical line node structure ~Py or a horizontal line of nodes in conjunction with significant gap modulation.Q 
. Since the question as to the nodal structure of the superconducting gap in Sr2Ru04 is still unresolved, we suggest 

■ here measurements of the anisotropy of the ultrasound attenuation a{Q) for rotating in-plane magnetic fields. Indeed, 
O ■ ultrasonic attenuation is another powerful tool for probing the anisotropic gap structure because a(0) is sensitive to 
O ■ the relative orientations of the sound direction and polarization, q and e, the field H(8), and the nodal directions of 

the gap. We shall show that a(G) exhibits fourfold symmetric variations in rotation angle & of the magnetic field while 
the variations of k[8] have twofold symmetry. It turns out that the field and temperature dependencies of a for vertical 
gap nodes are quite different for longitudinal and transverse sound modes while a for horizontal nodes is essentially 
. the same for longitudinal and transverse sound modes. Thus observation of the field and temperature dependence of 
" " ' a for longitudinal and transverse sound modes should yield important information on the nodal structure of the gap 
in Sr2Ru04. 

The ultrasonic attenuation a in the vortex state neaTrJIc2 for type-II s-wave superconductors has been calculated 
previously by S|Charnbergl3 and by Klimesh and Pesch.til In Ref. |l^ the Green's function G of Brandt, Pesch, and 
Tewordt (BPT)Ej has been employed which was derived from the spatial Gorkov integral equation for G(r, r', w) with 
kernel proportional to the " potential" 
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V^(ri,r2) = A{ri)A*{r2)exp[-2ie A ■ ds 
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where A(r) is the Abrikosov vortex-lattice order parameter and A is the vector potential of the magnetic field. 
This integral equation has been solved by expanding all functions in Fourier series k with respect to center of 
mass coordinates and in Fourier integrals p with respect to relative (difference) coordinates. In calculating physical 
quantities near Hc2 it often suffices to consider only the k = Fourier component G^o, lj). The corresponding 
anomalous Green's function F has been derived in the context of NMR theory by Pcsch.E^ The effective self energy 
part occuring in the denominators of these Green's functions is proportional to the spatial average of |A(r)p denoted 
by = A^pg(r) (1 — H/Hc2)- Furthermore it depends decisively on the quantity A/tisin6' where A = l/\/2eH 
, V is the Fermi velocity in the plane perpendiculat to H, and 9 is the angle between the quasiparticle momentum 
p and the field H. For 6^0 the Green's functions G and F tend to the ordinary Green's functions for an s-wave 
superconductor with gap A. For finite 9, the self energy contains both the effects of the Doppler shift and the Andreev 
scattering by the potential V^(p') = A^A^(5(p^)ea;p(— A^p'^)) for momenta p' in the plane perpendicular to H. This 
potential is multiplied by the hole-propagator G%p — p', — w) and integrated over all p'. The theory of a based on 
the Kubo formula and the BPT Green's functionll3 leads in the extreme clean limit to unphysical results in the final 
expressions. However, for small enough mean free paths £ this method is capable of describing experiments near Hc2 
sufficiently well. til 

Another approximation scheme for a, which is based on the Eilenberger equations and the Larkin-Ovchinnikov 
equations for the correlation functions via linear response theory,E3 can be carried out in the limit ^ — > oo. In this 
approximation scheme only the k = coefficient (spatial average) of the solution is considered. These authors state 
that this approximation is most questionable for quasiparticle directions p || H, and therefore they have concentrated 
on the special case of longitudinal ultrasonic attenuation where the wavevector q is parallel to H. Then the main 
contributions to the correlation function arise from directions p perpendicular to q, and thus perpendicular to H, for 
which their approximation is best. 

Since our main aim is to calculate a for unconventional superconducting gaps with nodes, we employ BPT Green's 
functions G and F where A^ is replaced by |A/(p)p with /(p) containing the p-dependence. This result has been 
derived from the original spatial integral equation for G'(r,r',cj) which contains in the kernal the non-local order 
parameters A(ri,r'i) and A*(r2,r'2) . The method consists in writing A(ri,r'i) = A(ri,ri — r'l) and A*(r2,r'2) = 
A*(r2, r2 — r'2) and introducing Fourier integrals with respect to the relative coordinates. In this way one can show 
that, to a good approximation {p' 1/A ^ pp), the A^ in the BPT Green's function is replaced by A^ |/(p)P- 

Another problem arises from the fact that we consider a for longitidunal and transverse sound waves with propa- 
gation vector q in the ab-plane of Sr2Ru04 , together with a field H in the plane which is rotated with respect to 
the a direction by the rotation angle 8. This means that the angle between q and H takes on all values, including 
q perpendicular to H so that the main contributions to a arise from directions of p near the direction of H. As has 
been pointed out above, for these directions the BPT Green's function is close to the ordinary Green's function which 
is physically quite plausible. Furthermore, for p ||H our correlation function, Eq. (H), tends to the pae derived in 
Ref. |l| for 7? = 0. Therefore we believe that, in contrast to the semi-classical approximation schemejlil our method 
can deal quite well with general directions of H with respect to q in the ab-plane. 

In Section II we present the general theory on the basis of the BPT Green's function. In Section III we discuss 
the ultrasonic attenuation results for different sound modes and several superconducting states with vertical and 
horizontal line nodes. The conclusions are given in Section IV. 

II. GENERAL THEORY OF ULTRASONIC ATTENUATION NEAR Hc2 

Our present theory for the ultrasonic attenuation a closely follows the method for calculating the thermal con- 
ductivity in magnetic fields near Hc2^ We start with the normal and anomalous BPT Green's functions G and F 
which corLtain both the effects of supercurrent ffow and scattering by the Abrikosov vortex lattice on the quasiparticle 
spectrum.Ej These Green's functions are employed in the expressions for the correlation functions for longitudinal 
and transverse sound propagation. Our method for evaluating the Kubo formula follows clpaely the method which 
has been used in the early theory of thermal conductivity by Ambegaokar and TewordtJlj First the integral of 
Re[G'(^, w) • G* (^ — V • q, — tJo) — F ■ F*] over the energy variable ^ is carried out. Here q and loq are the wave vector 
and frequency of the sound wave and v is the Fermi velocity. This yields the residue of this expression at the pole 
f = Co of the BPT Green's function G(C,a;) given by Eqs. (^) and (^. The most important term in the resulting 
expression is the quantity Im which yields the scattering rates due to impurity scattering and Andreev scattering 
by the vortices (see Eq. (^). We have obtained the Andreev scattering rate by calculating Im^o from Eqs. (||) and 
(H) for given w/A, AA/t), and 7 = (no impurity scattering) as a function of the angle 6 = (p — Q between the 
quasiparticle direction p and the field H. Then we find that for w > A (extended states) Im^o is zero in a range of 
angles above 9 = {p \\ H) which increases with lu, and becomes finite of order (AA/w)^ in a broad range of angles 
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up to § = Tr/2 (p _L H). (See Fig. (6)) It should be pointed out that the denominator of the second faiator in Eq. ^ 
(without the absolute square) occurs also in the expression for the density of states in the vortex stateJlj. In the limit 
9 = <j>—Q-^0{p II H) the expression for / in Eq. (^ (for q£ 0) tends to the well-known expression for the inverse 
relaxation rate Im^o Im(a)^ — |Ap)-^/^ (where uj — uj + i^) times the coherence factor 1 -I- (|(Dp — |Ap)/|(D^ — |Ap| 
which was first derived in Appendix C of Ref. In the hydrodynamic regime, ojqt <C 1, we obtain, by including 
vertex corrections in analogy to Ref. the following expression for the ratio of the ultrasound attenuation in the 
superconducting state, a^, to that in the normal state a„: 
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[l-^[AA/^sin(0-9)]^|/Pm^o)P] 
[(Imeo/r)2 + (g^)2(^.p)2] |l + 2[AA/i;sin(0-e)]2|/|2[l + z7^zoW'(^o)]P ' ^ ^ 

00 = (w + «7 + eo)[A/wsin(0 - 6)] , A = {2eH)-^'^ , 

7-r/.g, g^N{uj,H)/No (4) 

Here tt^^ and ir^y are the weight factors for the LlOO (TllO) and TlOO (LI 10) sound modes 00 where L — longitudinal, 
T = transverse, and q || [100] or [110] . We use the following notation: F is the normal-state impurity scattering rate, 
V is the in-plane Fermi velocity, 8 — Z(H, a) is the direction of H in the ab plane, A^ is the spatial average of the 
absolute square of the order parameter for the Abrikosov vortex lattice, and |/((?!>)P is the normalized absoLule square 
of the gap function. The quantity zq, and thus the pole of G, is given by the transcendental equationiEj'Ej 

zo = 2(w + ij)[A/vsm{(j) - 9)] + iy/^[AA/v sin(0 - 6)]^ \f\^ w{zo) . (5) 

Here, w{z) — exp{— z'^)er f c{—iz) . Note that Eqs. (^) and (||) differ from the corresponding equations in Ref. |l^ 
in that the sm.9 in the original quantity AA/t;sin0 is replaced by sin(0 — Q) for p lying in the direction </> and H 
lying in the direction G in the ab-plane, and A^ is replaced by A^ [/((/>) ^ for a gap with vertical lines of nodes. For 
a general state with /(p) one has to replace svnO with the expression given in Eq. (7) of Ref. |l^ and carry out the 
double integral over the polar and azimuthal angles and (/». In the following calculations we will neglect the term 
proportional to {qtf in the first denominator of Eq. (^) assuming the long wavelength limit g€ <C 1. The dependencies 
of (AA/w)^ and (A/w)^ on H/Hc2 are presented in Ref. |l^. The expression for ultrasonic attenuation in Eq. (^ is 
similar to the expression for k in Ref. |l^ apart from the missing factor w^, the weight factor 7r|^- instead of vf^ and 
the vertex corrections [second term in the curly brackets in Eq. (^]. For the vertex corrections we have used the 
expressions in Ref. |l^ for the unitary limit 6n ~ 7r/2 of the phase shift for impurity scattering. In the lirnit cj ^ , 
an analytical expression for the solution zq — i^o of Eq. (||) is obtained which yields, in close analogy to k Jij a much 
simpler expression for as /an in the limit T — > 0. 
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III. RESULTS AND DISCUSSION OF ULTRASONIC ATTENUATION FOR DIFFERENT SOUND 
MODES AND SEVERAL SUPERCONDUCTING STATES WITH VERTICAL AND HORIZONTAL LINE 

NODES 

First we consider the /-wave pairing state d = A z [px + 'iPy){Px ~ Py) for which |/|^ = cos^(20) has four vertical 
line nodes at </> = 7r/4, 37r/4, • • • on the cylindrical Fermi surface. In Fig. 1 we show our results for as /an versus 
H/Hc2 for the LlOO (9 = 0) and TlOO (6 = 7r/4) modes at T = and impurity scattering rate F/Aq 0.1. Aq 
is the BCS gap parameter, and we will always use the value (3a = 1-2 for the Abrikosov parameter. It is seen that 
as /an for the LlOO mode exhibits a strong upward curvature near Hc2 while it is a more linear function of H for the 
TlOO mode. For increasing F/Aq the upward curvatures decrease. 
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The solid curves for as/ an in Fig. 1 are calculated by neglecting the vertex corrections in Eq. (|^), while the dashed 
curves include the vertex corrections in the unitary impurity scattering limit. One recognizes that the effect of the 
vertex corrections in the unitary limit is rather small. This is also true for impurity scattering in the Born limit where 
in the denominator of Eq. (|l|) is replaced by one. It is interesting that, for field direction angles 8 = 7r/4 and 
= 0, the vertex corrections for the LlOO and TlOO modes, respectively, vanish. In the following calculations and 
figures we shall neglect the vertex corrections. 

In Fig. 2 we show the dependence of as/an on the in-plane field direction O at T = for fixed field strength, 
AA/u = 0.2, and impurity scattering rate F/Aq = 0.1. One sees that q;(0) has variations in of fourfold symmetry 
where the minima and maxima occur at = and 7r/4 for the LlOO mode. For the TlOO mode the maxima and 
minima are reversed. The maxima and minima can be explained as an effect of the 4> dependence of the weighting 
factors T:fj in Eq. and the density of states since the expression | • • • in the denominator of I{4>) in Eq. (j^) 
is proportional to | A'^((^, 0)p_ The density of states has a minimum (maximum) for quasiparticles traveling parallel 
(perpendicular) to the field.Ej The minima for LlOO and TlOO occur at = and 7r/4, respectively, because the 
corresponding weight factors tt^^. and ir'^y have maxima at = and 7r/4. For higher fields (e.g., AA/w — 0.1) a 
small local maximum and two neighboring minima occur around = 7r/4 for the TlOO mode due to the node of the 
gap at = 7r/4. 

The sound attenuation for the other /-wave pairing state, proportional to PxPy with |/p = sin^(20) , is obtained 
from the previous expression for |/p = cos^(20) by a simple transformation, 0' = — 7r/4, in Eq. (^. This exchanges 
the results for LlOO and TlOO which become functions of the new field rotation angle 0' = — 7r/4. 

We consider now the spin-triplet pairing state with horizontal line nodes □ where the squared gap amplitude is 
proportional to |/p = cos^ (cpz) = cos^[x] . Then, in addition to the cf) integration in Eq. (|l|), one also has to do the 
integration over x from — tt to -I-tt. In Fig. 1 we have included our result for as/an versus H / Hc2 for the mode LlOO 
at r = and F/Aq = 0.1. A strong upward curvature occurs near Hc2 which is similar to that of the LlOO mode for 
vertical nodes, but is quite distinct from the almost linear dependence of as /an on H/Hca for the TlOO mode in the 
case of vertical nodes. In Fig. 2 we have plotted our results for the in-plane field variation a{Q) for the LlOO and 
TlOO modes. Note that the function a{Q) for TlOO is obtained by shifting the function for LlOO by = 7r/4 along 
the axis. This can be seen by a variable transformation 0' = — 7r/4 in Eqs. (Q)"(^) making use of the fact that 
l/P = cos^(x) does not depend on (f). Comparison of the curves in Fig. 2 shows that the amplitudes of the variations 
a(0) for vertical and horizontal nodes are about the same (here about 10%). However, there is a marked difference 
in the form of the variations for the LlOO and TlOO modes for vertical nodes because, for the TlOO mode, the gap 
node manifests itself in the structure around 7r/4. In contrast to this, the form of the variations for horizontal nodes 
is the same apart from the shift by 7r/4 along the axis. 

We turn now to the lo dependence of the integrals in Eq. (|l|) which determine the T dependence of as/ a^ ■ It 
is now necessary to solve Eq. (||) for zq as a function of = w/A , and then to integrate the resulting expression 
for /((/)) in Eq. (^ (multipfied by irfjicl)), etc.) over cj) (and x for horizontal nodes). We denote by as{il.)/an the 
resulting integrand of the uj integral [without the factor (l/2r)sech^(Li;/2r)]. In Fig. 3(a) we plot as{^)/an vs 
for l/p = cos^(20) for large field magnitude [AA/v = 0.2 {H/Hc2 = 0.78)] and field directions = and n/A for 
both the LlOO and TlOO modes. It is interesting that the curves for = and 7r/4 cross at finite frequencies 
which means that the minima and maxima in a(0) are interchanged. Furthermore, the Q dependence for TlOO at 
= 7r/4 is almost linear while, for LlOO and = 0, it is quadratic. In Fig. 3(b) we show as{il)/an for the state 
|/|2 = cos2(cpz) . 

The temperature dependence of as /an is obtained according to Eq. (^ by integrating the expression 
(l/2T)sech2(w/2T) a,,(f^)/a„ over Lu . For the variable transformation from to, to 17 = uj/A we make use of the 
relations A/Aq = [1 - [H / Hc2)Y/'^ / y/JT^ and H/Hc2 = [1 + &l3A{AK/vf]-^ M In Fig. 4 we show our results for 
as /an vs T/Tc for vertical gap nodes at low and high field strength {AK/v = 0.6 and 0.2, or H/Hc2 = 0.28 and 0.78) 
and field directions = and = 7r/4 for the LlOO and TlOO modes. One recognizes that, for low T, the attenuation 
of the LlOO mode is approximately a quadratic function of T, while it is nearly linear for the TlOO mode at = 7r/4. 
These different low T power laws for the LlOO and TlOO modes agree roughly with the results of Ref. |lO| for zero field. 
It is interesting that the quadratic and linear T dependencies of as/ an for the LlOO and TlOO modes correspond to 
the quadratic and linear dependencies of as/a„ on H/ Hc2 for these modes (see Fig. 1). In Fig. 4 we also show as /an 
vs T/Tc for the LlOO mode for horizontal gap nodes and field directions = and 7r/4 . We recall that the results 
for LlOO for = (7r/4) are identical to those for TlOO for = 7r/4 (0) . One sees that the dependence on T is nearly 
quadratic. 
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IV. CONCLUSIONS 



In summary, we have calculated the ultrasonic attenuation a in layered Sr2Ru04 in the presence of in-plane magnetic 
fields for spin-triplet superconducting states with vertical or horizontal line nodes in the gap. This theory takes into 
account the effects of supercurrent flow and Andreev scattering (see Fig. (6)) by the Abrikosov vortex lattice near 
Hc2- For rotating in-plane field H{(d) the attenuation a{Q) exhibits variations of fourfold symmetry in the rotation 
angle Q. At T ~ the minima occur at = and 7r/2 for the LlOO sound mode, and at 8 = tt/A and 37r/4 for 
the TlOO mode. The amplitudes of the variations are about about the same (« 10%) for both horizontal and vertical 
nodes. However, in the case of vertical nodes, a{Q) for the TlOO mode shows characteristic structure in the directions 
of the nodes. For horizontal nodes the variations a{Q) for the LlOO and TlOO modes are the same if O is shifted by 
7r/4. The distinction between vertical and horizontal nodes also manifests itself in the different field and temperature 
dependencies of a. For vertical nodes the ratio as /an exhibits a strong upward curvature as a function of H/Hc2 near 
Hc2 for the LlOO mode while it is more linear for the TlOO mode. In the case of horizontal nodes, the attenuation for 
both the LlOO and TlOO modes shows a strong upward curvature near Hc2- The temperature dependence of the sound 
attenuation is determined by the frequency dependence of the integrand a(r2), il = ^/ in the expression for as/an- 
It is interesting that the curves for the functions a{Q) for field directions = and 7r/4 cross each other twice for 
increasing f2 indicating that the maxima and minima are interchanged twice. This causes corresponding crossings in 
the T dependencies of a for the two field directions. For vertical gap nodes, at low T, as/ an exhibits approximately 
a power law for the LlOO mode and a roughly linear T dependence for the TlOO mode. For horizontal gap nodes 
the functions as/an vs T/Tc for the LlOO and TlOO modes are identical for field direction differing by 7r/4. For low 
and high fields as /an follows roughly a power law. 

The ultrasonic, attenuation in Sr2Ru04 in zero field has been measured for the four in-plane modes LlOO, TllO, 
TlOO, and LllO.Q The attenuation follows a low temperature power law T^-^ for the LlOO (TllO) mode and r^ '* for 
the TlOO mode. In Ref. ^ a power law below T^^p and a dependence above T^^p has been measured for tlis 
TllO mode. Calculations based on the assumption of a circular cylindrical Fermi surface and vertical gap nodes t2l 
yield a linear T dependence of as/a„ for the TlOO (LllO) mode which disagrees with the measured T^ ** power law 
and, for the LlOO mode, a power law cloge to the measured one. The fact that the weakest T dependence occurs for 
the TlOO mode is taken as an indication Q that the gap nodes point in the [110] direction because this yields an excess 
of quasiparticles with wavevectors in this direction. 

Our results for the ultrasound attenuation for vertical nodes in magnetic fields up to Hc2 also yield the weakest 
(linear) H and T dependence for the TlOO mode because the corresponding weight factor has a maximum in the 
direction of a vertical gap node. Our calculations yield a stronger (quadratic) H and T dependence for the LlOO mode 
because the corresponding weight factor has a maximum in the direction of a vertical gap antinode. This is in strong 
contrast to the case of horizontal line nodes where we obtain a roughly quadratic H and T dependence for both the 
LlOO and TlOO modes, which are actually identical for in-plane field directions differing by 7r/4. 

At this point we would like to mention the modifications required if one takes the model of horizontal line nodes 
suggested in Ref. 0. This model can be simulated by multiplying the cos{cpz) used in this paper (see Ref. ^) by a gap 
with fourfold angular modulation, ( 1 -I- Acos4(/)), < |A| < 1, in the plane. This has the effect that the dashed curves 
in Fig. 2 for the variations a(9) of the LlOO and TlOO modes drift apart with a shift (and thus the anisotropy) that 
increases with increasing amplitude A of the fourfold gap modulation. See Fig. 5a for the case A = 1/2. Corresponding 
splittings occur in the attenuation curves a versus the field (Fig. 1) and temperature (Fig. 4) for the LlOO (TllO) 
and TlOO modes which are identical in our simple model of horizontal line nodes (see Fig. 5b). In fact, the H and 
T dependencies of a for the TlOO mode become weaker than those of the LlOO mode for increasing amplitude A 
of the gap modulation. We remark that such a fourfold angular modulation occurs in the gaps of the passive a— 
and /3— bands where horizontal line nodes are induced by the nodeless p-wave gap in the active band by interband 
proximity effect .tj We have fitted this gap modulation in the passive bands with our gap model discussed above and 
estimate a cather small value for A of about 0.1. The p-wave gap without line nodes is approximated in this model with 
A = —0.69. El For this p-wave gap the attenuation curves for the TlOO and LlOO modes in Figs. 1 and 2 (solid curves) 
are, roughly speaking, interchanged, which is plausible because the gap minima occur in the directions = 0, 7r/2, . . . 
. This has the effect that the H and T dependencies of a fot the LlOO mode become weaker than those of the TlOO 
mode. This appears to be in disagreement with experiment .0 

It should be pointed out that our theory is based on a number of simplifying assumptions. First, we have considered 
only the lowest order Fourier coefficient of the Green's functions with respect to the reciprocal Abrikosov vortex lattice. 
As the field decreases, the higher Fourier coefficients become more and more important and, near iJci, the Doppler 
shift method for a single vortex becomes more appropriate. Second, we have assumed a cylindrical single sheet Fernu 
surface, whereas the huge anisotropy observed between the TllO and TlOO normal state attenuation coefficients u 
indicates that the anisotropy follows from the nature of the three orbitals forming the three sheets of the Fermi 
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surface in Sr2Ru04.B3 These authors showed that in the case of Sr2Ru04 a simple viscosity tensor describing the 
electron-phonon interaction may not be sufficient for ultrasound attenuation. If this is the case our predictions could 
be changed. Third, we have assumed that the attenuation coefhcients normalized to their normal state values are 
insensitive to the observed normal state anisotropics. Here we have taken a temperature independent scattering rate 
r giving rise to aJjcmperature independent normal state attenuation which agrees roughly with the data for the LI 00 
and TlOO modes.!] Thus it cannot be excluded that the effect of the normal state anisotropy exceeds all anisotropics 
in the superconducting state calculated in this work. 

In spite of the above uncertainties we believe that measurements of the field dependencies (in addition to the 
temperature dependencies) of the ultrasonic attenuation for different in-plane sound modes can still provide useful 
information about the nodal structure of the energy gap in Sr2Ru04. 
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FIGURE CAPTIONS 

1. Ultrasonic attenuation ois/an vs applied field H/Hc2 aX T — Q for longitudinal and transverse in-plane sound 
modes LlOO (8 = 0; 9 = Z(H, a)) and TlOO (9 = 7r/4) for impurity scattering rate F/Aq — 0.1. Vertical gap nodes: 
solid and dashed curves without and with vertex corrections in the unitary impurity scattering limit, respectively. 
Horizontal gap nodes: dash-dot curve. 

2. asjoLn vs 9 = Z(H,a) for in-plane field rotation at T = for sound modes LlOO and TlOO and gap parameter 
AA/u = 0.2 {H/Hc2 — 0.78) and F/Aq = 0.1. Solid curves: vertical gap nodes. Dashes curves: horizontal gap nodes. 

3. Integrand Q;s(il)/a„ in Eq. (|l|) for Q!s/a„ (without the factor (l/2r)sech^(a;/2r)) vs reduced frequency = w/A 
for gap parameter AA/v = 0.2 {H/Hc2 = 0.78) and field directions 9 = (soHd curves) and 9 = 7r/4 (dashed curves), 
and F/Aq = 0.1. (a) LlOO and TlOO for vertical gap nodes; (b) LlOO for horizontal gap nodes. 

4. as/ttn vs T/Tc for vLlOO, vTlOO, and hLlOO (v = vertical, h = horizontal gap nodes), for field directions 
9 = (solid curves) and 9 = 7r/4 (dashed curves), (a) AA/v = 0.6 (-ff/iJc2 = 0.28), F/Ao = 0.2; (b) AA/v = 
0.2 {H/Hc2 = 0.78), F/Ao = 0.1. 

5. as /an for a gap with horizontal line nodes and fourfold angular modulation in the plane (A = 1/2) for the LlOO 
(sohd curves) and TlOO (dashed curves) sound modes, (a) as/a„ vs 9 = Z(H, a) for H/Hc2 = 0.78; (b) as/ an vs 
i//iJc2 at T = and for F/Ao = 0.1. 

6. Andreev scattering given by Im^o/^o (for 7 — 0), versus G = Z(p,H) for fixed A = AA/v « [{Hc2/H) — 
IY^'^/^/Wa and VL = uj/A. (a) A = 0.2 and f2 =0.0, 1.0, 1.2, and 1.5 (from top to bottom), (b) n = 1.0 and A = 
0.6, 0.2, and 0.1 (from top to bottom), (c) Vl = 1.2 and A = 0.6, 0.2, and 0.1 (from top to bottom). 
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